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Abstract 


A supercell storm is a tall, rotating thunderstorm that can generate 
hail and tornadoes. Two models exist for the development of the storm ’s 
rotation or mesocyclone — the conventional splitting-storm model, and 
the more recent pulsing inertial oscillation (PIO) model, in which a 
nonlinear pulse represents the supercell. Although data support both 
models and both could operate in the same supercell, neither model has 
satisfactorily explained the tornadic cyclone. A tornadic cyclone is an 
elevated vorticity concentration of Rossby number Ro » 1 000 that 
develops within the contracting mesocyclone shortly before a major 
tornado appears at the surface. We now show that if the internal 
temperature excess due to latent energy release is limited to the realistic 
range of -12 K to +12 K, the PIO model can stall part way through the 
pulse in a state of contraction and spin-up. Should this happen, the 
stalled-PIO model can evolve into a tornadic cyclone with a central 
pressure deficit that exceeds 40 mb, which is greater than the largest 
measured value. This simulation uses data from a major tornadic 
supercell that occurred over Oklahoma City, Oklahoma, USA, on May 3, 

1999. The stalled-PIO mechanism also provides a strategy for human 
intervention to retard or reverse the development of a tornadic cyclone 
and its pendant tornado. 

1. Introduction 

Supercell storms are tall rotating thunderstorms that can generate hail and tornadoes. Doppler radar 
images of a hailstorm (Miller, Tuttle, and Knight 1988) show a complex structure that includes an updraft 
that peaks at midheight, two downdrafts, and a mesocyclone. These features evolve in an environmental 
flow whose velocity increases and rotates with height. The environmental flow includes a boundary layer 
at the surface, a highly sheared layer near the storm top, and a midtropospheric layer in which the envi- 
ronmental flow is blocked and must circumnavigate the storm (Brandes 1981). Our modeling study, 
which is still in a rudimentary state, is limited to the midtropospheric layer and assumes that compatible 
flows exist in the boundary and upper layers. 

Latent energy release drives the vertical flow; however, the source of the rotational flow or mesocy- 
clone is controversial. Two models exist for the generation of the mesocyclone. The first is the conven- 
tional splitting-storm model, as reviewed by Doswell (2001), Church et al. (1993), Ray (1986), Klemp 
(1987), Kessler (1986), Bengtsson and Lighthill (1982), Schlesinger (1980), and others. The second is 
the pulsing inertial oscillation (PIO) model (Costen and Miller 1998 and Costen and Stock 1992). Data 
exist that support both the conventional and the PIO models, and both mechanisms could operate in the 
same supercell. 

In the conventional model, a buoyant updraft tilts and stretches horizontal vortex lines in the vertically 
sheared environmental flow to produce two counterrotating supercells that separate. Such splitting super- 
cell pairs are observed. In fact, we will use data from the cyclonic member of a splitting pair that oc- 
curred on May 3, 1999 over Oklahoma City, Oklahoma, USA. However, Davies-Jones (1986) states that 
observed supercells are predominantly cyclonic. In his modeling studies, Klemp (1987) explains this 
apparent discrepancy by showing that environmental conditions documented by Maddox (1976) will often 
weaken or suppress the anticyclonic member of the pair. 



Data supplied by Davies-Jones (1986) show that the generation of a major tornado requires a transition 
wherein the mesocyclone, or a portion thereof, contracts and spins up to a Rossby number Ro ® 1000, 
where Ro is the ratio of the vertical spin rate of the cyclone core to the local vertical spin rate of the Earth. 
This contraction first occurs in the midtropospheric layer, as observed by Doppler radar data presented by 
Davies-Jones (1986), Vasiloff (1993), and Burgess and Magsig (2000). The contracted mesocyclone is 
termed a tornadic cyclone because its development aloft precedes the onset of the strongest tornadoes at 
the surface. 

Although the conventional model is well developed, it has not satisfactorily simulated a tornadic 
cyclone. Modeling studies by Wicker and Wilhelmson (1995) and Wilhelmson and Wicker (2001) pro- 
duced vorticity concentrations that are strong enough for a tornadic cyclone; however, these concentra- 
tions originated at the surface instead of aloft — in conflict with the data. Modeling studies of the flow in 
a rotating cup by Trapp and Davies-Jones (1997) showed that the height where tornadic vorticity first 
occurs depends on the vertical distribution of buoyancy; however, it is not clear that their results apply to 
a free mesocyclone in the troposphere. 

The PIO model is in its infancy and does not yet resolve individual updrafts and downdrafts but 
instead treats the average vertical flow in the core — represented by a uniform time-dependent vertical 
flow. The model is essentially a Rankine vortex that has been generalized to include the Coriolis force 
plus this uniform vertical flow within its circular cylindrical core. Ferrel (1889) attributed supercell rota- 
tion to the Coriolis force. Scaling studies by Morton (1966), Wallace and Hobbs (1977), Davies-Jones 
(1986), Holton (1992), and others have since concluded that the Coriolis force is too weak to spin up a 
supercell storm in the observed time. However, scaling studies rely upon the mechanism, and none of 
these studies considered the PIO mechanism. 

In the compressible continuity and momentum equations used in the PIO model, the vertical tempera- 
ture lapse of the middle troposphere is neglected and the density is taken to decrease exponentially with 
height. For our purposes, this approximation is better than the constant density approximation used in the 
analytic studies of supercell storms by Rotunno and Klemp (1982) and Filly (1986). However, we will 
retain the actual temperature lapse in the energy equation and for stability considerations of the midtropo- 
spheric layer. Our practice of neglecting terms where they are presumed to be unimportant and retaining 
them where they are known to be important has precedent elsewhere in Atmospheric Science, for exam- 
ple, in the Boussinesq approximation (Menke and Abbott 1990), which we do not use. 

In our postulated environment, rising parcels of air expand (dilatate) and descending parcels contract. 
The key approximation used in the PIO and stalled-PIO models is that the horizontal divergence in the 
core results predominantly from this three-dimensional dilatation; that is, the horizontal divergence is 
approximately equal to the speed of the vertical flow divided by the scale height. We shall refer to this 
approximation as the “dilatation-horizontal divergence” approximation or simply the “DHD” approxima- 
tion. It follows that when the vertical flow in the core is downward, the cylindrical core contracts. This 
contraction causes the core to spin up, partly by increasing the concentration of existing vertical vorticity 
and partly by the action of the Coriolis force on the convergent flow. Fikewise, upward flow causes the 
core to expand and spin down. This feature of spin-up occurring in a downdraft is in contrast to the con- 
ventional splitting- storm model, where spin-up occurs in an updraft. 

Various arguments can be advanced to support the DHD approximation, as in the statement that fol- 
lows equation (2.4); however, it seems clearest simply to test the accuracy of the DHD approximation 
during the model runs, and this we will do. The DHD approximation neglects certain terms in the conti- 
nuity equation. We will compute both the retained terms and the neglected terms and will show by plot- 
ting their ratio that the DHD approximation is reasonably accurate during our model runs. 
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In the PIO model, the generalized Rankine vortex contracts, expands, and translates in the midtropo- 
spheric layer in accordance with a nonlinear inertial oscillation, whose period is typically about 19 hr. 
The pulse phase of this oscillation, which coincides with the contracted state of the vortex, is identified 
with a single cyclonic or anticyclonic supercell. 

According to data presented by Costen and Miller (1998), only the pulse phase of the PIO model is 
actually observed in the troposphere; that is, the oscillation starts at the beginning of a pulse, proceeds 
through the pulse (supercell storm), and terminates shortly afterwards. Therefore, we will limit our pre- 
sent simulations to the pulse phase. For simulations of the complete period, see Costen and Miller (1998) 
and Costen and Stock (1992). 

Because of the DHD approximation and the prescribed vertical dependence of density, the continuity 
and momentum equations are sufficient to obtain the PIO and stalled-PIO solutions for the vorticity, hori- 
zontal divergence, velocity, pressure, core radius, core buoyancy, mesocyclonic circulation, and cloud 
base mass influx. The equation of state for a diy, perfect gas then determines the temperature, which, as 
mentioned, turns out to be isothermal, except for a uniform temperature excess or deficit in the core. The 
energy equation serves only to determine the thermal input power or latent energy release that is required 
to support the PIO and stalled-PIO solutions. 

Although the PIO is driven by latent energy release, the oscillation is controlled by horizontal inertial 
flow (Haltiner and Martin 1957 and Holton 1992) that is organized into a radial oscillation. When viewed 
from above (fig. 1(a)), every parcel on the periphery of the contracting core traverses an anticyclonic 
circular arc such that the parcel’s centrifugal force balances the Coriolis force and the horizontal pressure 
gradient is zero. Costen and Stock (1992) show that the same is true for every parcel inside the core. 
Consequently, the PIO model has no central pressure deficit (prior to a stall). During contraction, each 
parcel in the core is compressed by descending in a uniform downdraft. When the core reaches its mini- 
mum radius and maximum spin-up, the downdraft ceases. The subsequent inertial trajectories (fig. 1(b)) 
cause the core to expand and spin down, which requires an upward flow in the core. 

Compared with the complicated structure of a supercell described earlier, the PIO model, in its present 
state of development, is a zeroth-order approximation. The vertical speed, horizontal divergence, vortic- 
ity, and buoyancy — all of which are uniform in the core of the PIO model — are zeroth-order approxima- 
tions to the instantaneous spatial average values of these fields in the core of an actual supercell. The 
vertical flow reversal (described in the previous paragraph) that occurs at the midpoint of the pulse in the 
PIO model therefore represents a reversal in the average vertical flow of an actual supercell. 

The PIO model determines the amount of core buoyancy required to sustain the oscillation. The aver- 
age downdraft that occurs during spin-up at the onset of a pulse requires a moderate amount of negative 
buoyancy that results from evaporative cooling of midtropospheric moisture. However, the average draft 
reversal that occurs at the midpoint of the pulse requires a large spike of positive buoyancy. This spike is 
generated by condensational heating of moist surface air that is lifted in a convectively unstable environ- 
ment by the gust front of the downdraft. 

Costen and Miller (1998) showed that the available buoyancy range limits the maximum Rossby num- 
ber that the PIO can achieve. We define the available buoyancy by a core temperature excess A T in the 
range (-12 K < AT < 12 K), as suggested by sounding data presented by Miller, Tuttle, and Knight 
(1988). The PIO solution then gives a maximum Rossby number of Ro ~ 100, which is an order of mag- 
nitude too low for a tornadic cyclone. However, we will now show that this same available buoyancy 
range can cause the PIO model to stall in a state of contraction and spin-up. Our objective is to model the 
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flow after such a stall has occurred and to demonstrate that it can produce a tornadic cyclone. Since the 
PIO and stalled-PIO models apply in the midtropospheric layer, the resultant tornadic cyclone develops 
aloft, in agreement with observation. 

A remarkable feature of the PIO model is that its track over the surface of the Earth during the pulse 
phase resembles the dogleg track of an actual supercell hailstorm, as shown by Costen and Miller (1998) 
and Costen and Stock (1992). After a stall has occurred, however, the track equations disappear from the 
model, so we will not attempt any track simulations here. The issue of the post-stall track must await 
further research. 

As mentioned, the PIO and stalled-PIO models will use data from the May 3, 1999 supercell storm 
that occurred over greater Oklahoma City, Oklahoma. This supercell (Supercell A) was the right-hand 
member of a splitting supercell pair. It generated a 1.6-km-wide tornado (Tornado A9) that killed 
36 people and injured 583. This tornado also inflicted property damage of 1 billion United States dollars 
(USD) along its 60-km track. Rated F5 (maximum intensity) on the Fujita scale (Fujita 1973), the tor- 
nado possibly reached wind speeds of 3 1 8 mph, as indicated by preliminary Doppler radar data. Burgess 
and Magsig (2000) present time-height data that confirm that an elevated tornadic cyclone developed in 
Supercell A about 5 min before the start of Tornado A9, that this tornadic cyclone persisted with some 
variation throughout the 80-min lifetime of Tornado A9, and that the tornadic cyclone dissipated about 
3 min before this tornado ended. 

Although we use data from Supercell A, we will not attempt a precise simulation of this storm. In- 
stead we will explore the capabilities of the PIO and stalled-PIO models within the framework of these 
data. While awaiting publication of a value for the mesocyclonic circulation of Supercell A, we have 
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used a typical value given by Davies-Jones (1986) of T max =5x10 m“ s . We will always initialize 
the core radius such that this value for circulation is achieved at the time of maximum contraction 
and spin-up. As mentioned, the model computes the thermal input power Q(t) from latent energy 
release and the mass influx at cloud base M(t). While also awaiting these values from Supercell A, 

13 

we have compared our computed values with measured reference values of <2 re f = 1.9x10 W and 
9 —1 

M re f = 1.2 x 10 kg s presented by Foote and Fankhauser (1973) for a different supercell. 

Sections 2 to 4 give a review of the PIO model and point up a naturally occurring condition that could 
cause the model to stall in a state of contraction and spin-up. Sections 5 and 6 give the mathematics of 
tornadic cyclone development after such a stall has occurred. In section 7, the stalled-PIO model is 
applied to Tornadic Cyclone A9, with some attention given to the possibility of human intervention to 
promote early termination. Concluding remarks are presented in section 8. 

2. Relevant Aspects of PIO Model 

2.1. Tangential Coordinate Frame 

The model uses the local tangential Cartesian frame shown in figure 2. The origin is fixed at mean sea 
level (MSL) on the Earth’s surface, and the x, y, and z axes point eastward, northward, and upward, 
respectively. (Cylindrical coordinates r, 0, and z are also used, where x = rcos0 and y = r sin 0.) 
Although the x and y axes do not curve with the Earth’s surface, Haltiner and Martin (1957) show that 
this frame is adequate for describing tropospheric flows within a horizontal radius of about 60 km. This 
radius is sufficient to contain the convective core of a supercell, but not all the outer flow. In this frame, 
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the fluid velocity is denoted by \ = (u,v,w) (m s -1 ), the vorticity by ct> = curlv = (|,r],5) (s-l), the 
divergence by D = divv = du/dx + dv/dy+dw/dz (s~l), and the angular velocity of the Earth by 
Q = (0,Q y ,Q_) (rads- 1 ). The components Q r and Q„, although slowly varying functions of y, are 
treated as constants, which is reasonable for the horizontal radius given previously. 

2.2. Generalized Rankine Vortex 

The generalized Rankine vortex shown in figure 2 is intended to represent the convective core of a 
supercell storm that upon further contraction and spin-up can become a tornadic cyclone. The core radius 
is a(t). The centerline is located at [x c {t), y c (t)\. The vertical vorticity l,(?) is uniform inside the core 
\r < n(t)] and zero in the outer region [ r > a(t)]. The midtropospheric layer, which is between the 
planetary boundary layer below and the highly sheared layer above, is given by ( b<z< h). The density 

p (kg m -3 ) is taken to decrease exponentially with height z. This stipulation neglects the temperature 
lapse between b and h. Rising/falling parcels of air dilatate and contribute to the divergence D(t), 
which is also taken to be uniform inside the core and zero in the outer region. The core fluid is uniformly 
buoyant, with a normalized density deficit B(t), which is called simply the buoyancy. Lateral entrain- 
ment, friction, and heat conduction are neglected, which is reasonable for a supercell storm in the mid- 
tropospheric layer. 

The density of the outer fluid at radius a is given by 

p°(a,z) = pJe _8Z (2.1) 

where the fields in the outer region are distinguished by the superscript 0 and Z = z-b. The coefficient 
p° h and the inverse density scale height e (in m-1) are constants. Inside the core, the density is given by 

p(z,0 = pJ[l-5(0]e- £Z (2.2) 

The buoyancy B(t)« 1, because it corresponds to a temperature excess or deficit of at most 12 K, as 
mentioned, and terms of 0(B ) will always be neglected. 

The continuity equation is given by 


— + v • Vp + pD = 0 (2.3) 

dt 

Substituting the inner density (eq. (2.2)) in equation (2.3) and dividing by p gives 


-B -Ew+ D = 0 (2.4) 

where an overhead dot denotes an ordinary time derivative. Retaining the effects of B(t) on the buoyant 
force but neglecting its effects on D(t), we have 

D(t) = £w(i) (2.5) 
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which is the key approximation used in both the PIO and stalled-PIO models, wherein D and ew will be 
used interchangeably. Equation (2.5) confirms that the vertical component of velocity w(t) is also 
spatially uniform inside the core so that D(t) is purely horizontal divergence given by 


_ . . du dv 
D(t) = — + — 
dx dy 


( 2 . 6 ) 


Hence, approximation (eq. (2.5)) states that the horizontal divergence inside the core results 
predominantly from the dilatation of ascending/descending air parcels. We call equation (2.5) the 
“dilatation-horizontal divergence” approximation, or simply the “DHD” approximation. To check the 
accuracy of this approximation, we will plot the ratio i? DHD approx neglected to retained terms in 

equation (2.3); for the PIO model this ratio is given from equation (2.4) by 


„no 

^DHD approx 


B 

D 


(2.7) 


Also included in the core is a horizontal flow with vertical wind shear that is represented by the constants 
U, and V, and the horizontal vorticity components t(/) and i](r). 

3. Inner Reduction 

3.1. Velocity Field 

The complete velocity field inside the core is given by 

u(x, y,z,t) = U b + r\(t)Z + ^ [D(t)X(x, t) - 2,(t)Y(y, t)] (3.1a) 

v(x, y,z, t) = V b - £(, t)Z + ^ [D(t)Y(y, t) + £(0 X(x, t)] (3.1b) 

w(x,y,z,t) = w(t) = (3.1c) 

e 

where X(x,t) = x-x c (t) and Y(y,t) = y - y c (t). The constants U f) and V t) represent the stationaiy part 
of the horizontal flow inside the core, and r|(t)Z and - ^(t)Z represent the inner vertical wind shear. 
The terms -C(t)Y/2 and 2(1) X/ 2 represent the inner flow of the Rankine vortex. The angular velocity of 
the core fluid 0 (in rad s _1 ) and the Rossby number Ro are given by 

• t t 

0 = — Ro = — (3.2) 

2 2Q z 

The terms D(t)X/2 and D(t)Y/2 represent the inner radial horizontal flow that results from the 
horizontal divergence. 
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3.2. Momentum Equation 


The inviscid, compressible momentum equation is given by 

— + coxv + V — + — + V<l> + 20 xv = 0 (3.3) 

dt 2 p 

where v 2 is v • v, p is the pressure (Pa), O is the geopotential gz (m- s~ 2 ), and g is the gravitational 
acceleration (9.81 m s~ 2 ). 

Solving equation (3.3) for V/; and substituting the inner velocity field (eqs. (3.1)), we find 


where 


dp 

dx 

dp 

dy 

dp 

dz 


P 

2 

P 

2 


- G\i’+ Duj 
(G., + (Xu + Dv) 

( D \ 

g + 2Q u 

P s 


Gj = DX - 1 Y-Dx c + ly c + 2i]Z + 



(3.4a) 

(3.4b) 

(3.4c) 

(3.5a) 


G,=DY+iX-lx-Dv-2tZ-2 — l (3.5b) 

e 

G 3 = ^ + 4Q _ (3.5c) 

3.3. Second-Order Partial Derivatives and Nonlinear Harmonic Equations 

In the midtropospheric layer, we assume that p and its first- and second-order partial derivatives are 
continuous functions of x,y, andz (except at radius a, where p must be continuous but its derivatives 
could be discontinuous). It follows for (r <a, b<z < h) that 

PI P2 PI P2, 

dp dp dp dp 

dx dy dy dx dx dz dz dx 

Substituting the cross derivatives of equations (3.4) in equations (3.6), dividing by p, and setting the 
coefficients of X, Y, and Z individually to zero in each equation, we obtain the following set of coupled 
ordinary differential equations (ODEs): 

i-i[il(S + 4Q z )-Z>i] (3.7a) 


P2 

d p 


P2 
a p 


dy dz dz dy 


(3.6) 
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f] = -^[^ + 4Q z ) + Z)ri] 


(3.7b) 


t = -Z)(^ + 2Q z ) (3.7c) 

J D = i^(C + 4Q_„)-D 2 ] (3.8) 

x c =U b + - Q v + ^ ~ 2 \^(-^ w+2Q z U b) + D (^ w - 2Q z V b)\ ( 3 - 9a ) 

8 L, + D" 

h- r i,*^^[i(-^* 2Q , v b)* D (-^* 2Q , u t)] < 3% > 

3.4. Pressure Field 


To obtain the inner pressure field p, we first substitute equations (3.7), (3.8), (3.9), and (2.2) into the 
pressure gradient (eqs. (3.4)), which gives 

^ = -p 9 e- £Z (l -B)Q - (3.10a) 

dx * e 

^ = pJe- EZ (l-£)Q v ^ (3.10b) 

dy - £ 

^ = -p 9 e- 8Z (l-5)|g + w-2Q j; U b + ^Z + ^(DX-IY) J (3.10c) 

These partial derivatives can now be integrated to obtain p as 

2q 

p(x,y,z,t) = P y £Z (l-B)Z 1 + - y. I\. + U b + r { Z + -(DX-l Y ) (3.11) 

e g g L 8 2 J 

This result confirms that the PIO model has no pressure deficit at its center, where X = Y = 0, which is 
consistent with inertial flow. According to Davies-Jones (1986), the pressure deficit in mesocyclones is 
typically about 5 mb, although it has been measured as high as 34 mb. This discrepancy indicates the 
occurrence of a partial transition from inertial flow to cyclostrophic flow, where the centrifugal force on 
parcels is balanced by the radial pressure gradient, as discussed in sections 5 and 6. 

Thus far, we have reduced the inner solution to a time-dependent set of ODEs, except that we have no 
equations for a(t ) or B(t). These quantities are determined by the jump conditions at the cylindrical 
interface between the inner and outer solutions, as shown in figure 2. Instead of attempting a nonsymmet- 
ric, translating outer solution, we will revert to an axisymmetric, nontranslating solution that is sufficient 
for obtaining a{t) and B(t). We assume that the results so obtained apply approximately to the 
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translating model. This assumption will fail if the inner flow of the translating model does not block the 
relative environmental wind, which will then erode the core. 

4. Axisymmetric, Nontranslating PIO Model 

4.1. Inner Solution 

For the PIO model to remain centered on the origin (i.e., for x =j =0), we must take 
U b = V b = | = r| = Q = 0. The velocity (eqs. ( 3 . 1 )) and pressure (eq. ( 3.11 )) then become, in cylindrical 
coordinates, 


, , D(t) 

v r (r,t)= r 

(4.1a) 

v 0 (r,O= 2 r 

(4.1b) 

v (t) = w(t) = 

E 

(4.1c) 

p(z,0 = pJe- EZ ^[l-f?(t)][l + ^ 

E L s 

(4.2) 

4.2. Outer Solution 


Because the external vorticity co° and divergence are zero, 

(for r > a) 

the outer velocity field is given by 

2 r 

(4.3a) 

v „ tiOa'-U) 

2 r 

(4.3b) 

II 

O 

II 

o 

(4.3c) 


The field v } , is continuous at the interface r = a(t), as required by the jump conditions for the conti- 
nuity and momentum equations. The field v 0 is also continuous, but w, p, D, and t, all have finite dis- 
continuities. The requirement that the interface moves with the fluid gives the desired equation for core 
radius 



(4.4) 
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The core buoyancy B(t) can be determined from the jump condition that the pressure is continuous at 
r = a. To determine the outer pressure field, we first solve the momentum equation (3.3) for the pressure 
gradient in cylindrical coordinates and then substitute the outer velocity field (eqs. (4.3)) and the tenden- 
cies (eqs. (3.7c), (3.8), and (4.4)) to obtain 


3 0 
dp 

dr 


gF P° 


2\ 


1 -- 



(4.5a) 


dp 

ae 


0 


(4.5b) 


. 0 
dp 

dz 



(4.5c) 


where F, the Froude number, is a constant of the motion given by 

F = ^l(^ +D 2\ 

4 g \ I 


(4.6) 


Its constancy for the PIO model is demonstrated by taking its time derivative and substituting the 
tendencies (eqs. (3.7c), (3.8), and (4.4)). 

As with the inner solution, we require 


.20 .20 

d p _d p 

dr dz dz dr 


(r > a, b<z< h) 


(4.7) 


Consequently, the outer density p u must satisfy 


dp° _ F ( 
dr er 


1 -- 


dp 

dz 


0 


(4.8) 


The solution of this equation that satisfies the boundary condition (eq. (2.1)) at r = a is given by 


p°(r,z,t) = pj 



F 

„ F 

( 2 ) 



exp 

-eZ + — 

1 


J) 


2 

r 2 

\ 1 y 



(4.9) 


Substituting this result in equations (4.5) and integrating for the outer pressure, we obtain the hydrostatic 
result 


P°(r,z,t ) 


^p°(r,z,0 = p°- 

£ 8 


/ \ 

F 


{ 2) 


a 

exp 

-eZ + — 

1-^ 


J, 

2 

i S') 



(4.10) 
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Both p°(r,z,t) and p °(r,z,t) have their maximum values at r-a and approach zero as r —* <», 
which is reasonable for our tangential coordinate frame. To conform with the notation for the stalled-PIO 
model that will be developed in sections 5 and 6, we associate the maximum external pressure p°(a,b) 
and density p °(a,b) at height b with the environmental pressure /;^ nv " and density p^ nv " at height b, 
which are constants. Then, by equation (4.10), we have 


0 _ envir _ e envir 

Pb~Vb ~~ P b 

o 


(4.11) 


and the expressions for the outer pressure (eq. (4.10)), inner pressure (eq. (4.2)), and inner density 
(eq. (2.2)) become 


P° (r, z, t) = - p° (r, z,t) = pf 
e 


(a) 

F 

„ F 

{ 2 1 



exp 

-eZ + — 

1 


J, 

2 

r 2 

) ' ( 



(4.12) 


M Z ,0 = 4 nv V £Z [1-5(0] 


1 + 


w{t) 


g 


(4.13) 


P {z,t) = -p e r^-B) e- eZ 
g 


(4.14) 


4.3. Buoyancy and Inner Pressure 

By equating the outer pressure (eq. (4.12)) and the inner pressure (eq. (4.13)) at r = a, we obtain the 
desired equation for buoyancy 


B(t) = 


Mt) 

g 


(4.15) 


This expression is similar to that given by Darkow (1986) for the vertical acceleration of a nonentraining 
buoyant parcel. The inner pressure is now independent of t and is given by 

p(z)-4™e-' Z (4.16) 


By equation (3.8), the buoyancy (eq. (4.15)) for the PIO can also be determined from the formula 


B no 


1 ' 

2 Eg ■ 


t(^ + 4Q z ) 



4.4. Supplementary Formulas 

The mesocyclonic circulation T is given by (in m 2 s _1 ) 

T = n a 2 t, 


(4.17) 


(4.18) 
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and the cloud base mass influx M by (in kg s -1 ) 


M = 2jt w Jq rp(r,b,t ) dr (4.19) 

Substitution of equation (4.14) gives for the PIO 

M n0 =-a 2 D(l-B)pl n ™ (4.20) 

g 

The temperature is determined from the equation of state for a dry perfect gas 


p = pRT (4.21) 

which we apply to both the inner and outer regions, where R, the gas constant, equals 287 J kg -1 K -1 . As 
mentioned, the outer temperature 7’° in the model is a constant 


T 


0 


g_ 

eR 


(4.22) 


In practice, this formula determines the inverse scale height e from a measurement of T ° , which we will 
take to be the environmental temperature at height b. The inner temperature T is given by 


r(0 = T°[l + 5(0] 


(4.23) 


Thus, the model core is also isothermal, but its temperature varies with the buoyancy. Solving equa- 
tion (4.23) for B , we confirm that 


T(t)-T° _ AT(t) 
1 ; T 0 


The energy equation for a dry perfect gas, as given by Holton ( 1 992), can be written 


<l = P c p 


(d J 


fd J 

[dt J 

T- 

— + v-V 
[di J 


P 


(4.24) 


(4.25) 


where q (in W rrr 3 ) is the diabatic heating rate, which is usually due to latent energy release, and c p 

equals 1004 J kg -1 K -1 . Substituting the inner pressure (eq. (4.16)) and temperature (eq. (4.23)) and 
integrating over the core volume, we obtain the total thermal input power Q (in W) required to support 
the PIO 


Q 


PIO 


envir 

~ Pb 


-cT°B + D 

MV,)' 

[l-e“ B M 

g P 

U * ). 



(4.26) 
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As explained, although we have neglected the temperature lapse dT°/dz in the continuity and momentum 
equations, we have retained it here to be more accurate when applying this formula to an actual supercell. 
The corresponding water vapor influx M y (in kg s _1 ) is given by 



where L , the latent heat of condensation, equals 2.5 xlO 6 Jkg -1 . 

4,5. Maximal PIO Plots 

Plots of A T, w, a, Ro, R DHD approx , T,M, and Q versus t are shown in figure 3 for a 2-hr interval 

centered on the pulse phase of the PIO. The corresponding inertial trajectories of parcels on the core 
periphery were previously shown in figure 1(b). Although equations (3.7c) and (3.8) have an analytic 
solution (Costen and Miller 1998), these two equations and equations (4.4) were integrated numerically 
by using a fourth-order Runge-Kutta routine. The relevant parameters, listed in table 1, are taken from 
data on Supercell A that was provided by Edwards and Thompson (2000) and Burgess and Magsig 
(2000). The initial values listed in table 2 were chosen so that the mesocyclonic circulation at maximum 

contraction and spin-up has the value T max =5x10 ms and the core temperature excess AT plotted 
in figure 3(a) just stays within the available range (-12 K < AT < 12 K), which corresponds to buoyancy 
B{t) in the range (-0.04158 < B(t) < 0.04158). 

As shown in figure 3(a), AT peaks at +12 K, which corresponds to the maximal PIO. The corre- 
sponding core buoyancy drives the vertical speed w shown in figure 3(b). By the DHD approximation 
(eq. (2.5)), the horizontal divergence D is proportional to w, so the core radius a varies as shown in 
figure 3(c). The resulting peak Rossby number, as shown in figure 3(d), is the largest the PIO mechanism 
can produce for the environment of Supercell A. As mentioned, this value Ro = 118.5 is an order 
of magnitude less than that required for Tornadic Cyclone A9. Figure 3(e) is a plot of the ratio 
^DHD approx’ as given by equation (2.7), and shows that this ratio is very small except for a brief interval 

in the middle of the pulse when it reaches 0.12 as both B and D pass through zero. We conclude that 
the DHD approximation was reasonably accurate during this PIO model run. 

The mesocyclonic circulation T, as shown in figure 3(f), starts out anticyclonic, achieves the targeted 

5 2—1 

cyclonic peak value of r max =5x10' m s , and eventually becomes anticyclonic again, in agreement 

with the inertial trajectories shown in figure 1(b). The cloud base mass influx M from equation (4.20) is 
shown in figure 3(g) and the thermal input power Q from (4.26) in figure 3(h). The perturbation in Q at 

r = 60 min results from the B term; otherwise, Q is negative during the downdraft phase and positive 
during the updraft phase because downdrafts generally require evaporative cooling and thus are exother- 
mic, while updrafts require condensational heating and are endothermic. 

A comparison of figures 3(a)-(d) confirms that spin-up occurs in a contracting cylindrical downdraft 
that is driven by moderate negative buoyancy. At the midpoint of the pulse, the contracting downdraft 
reverses and becomes an expanding updraft that causes the mesocyclone to spin down. This rapid draft 
reversal requires a large spike of positive buoyancy. 
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4.6. Condition for PIO To Stall 


Figure 3(a) depicts an 8 to 1 asymmetiy between the +12 K and -1.5 K excursions of the core tem- 
perature excess AT. Since the negative value of AT occurs first, we are free to initialize the PIO model, 
as shown in table 3, so that AT decreases to a minimum value of -12 K. The subsequent maximum 
value of +96 K is far outside the given available range of (-12 K < AT <12 K). We infer that upon 
reaching AT = +12 K, the PIO model would stall, as shown in figure 4(a). Our task now is to develop a 
model that applies after such a stall has occurred. 

The physical mechanism shown in figures 1 and 3 gives insights about what should happen after a 
stall. The truncated positive buoyancy shown in figure 4(a) is insufficient to cause the rapid draft reversal 
required for the core to expand in accordance with the inertial trajectories shown in figure 1(b). Conse- 
quently, the parcels would become trapped in orbits about the contracted cyclone center. Since the 
Coriolis force could no longer balance the centrifugal force, depressions in the pressure and density would 
begin to develop at the center, and the flow would undergo a transition from inertial to cyclostrophic. 

5. Stalled-PIO Solution 


To allow for the development of such depressions in pressure and density after the stall, we must gen- 
eralize the outer and inner densities as follows. At the interface r = a, we now take 

p °(a,z,t) = o° h (t)e- EZ (5.1) 

p(a,z,t) = o° h (t)[l-B(t)]e- £Z (5.2) 

and for the inner density 

p(x, y,z,t) = o(X,Y,Z,t)[\-B(t)] e" EZ (5.3) 


where B(t ) is now a given function of 1 that is confined to the available range (-0.04158 < 
B(t) < 0.04158), which corresponds to (-12K < AT(t) < 12 K). Although the inner analysis will use 
o(X,Y,Z,t) to establish certain generalities, the model runs will revert to a simplified axisymmetric 
o(r,t). Substituting this simplified inner density into the continuity equation (2.3) gives 


a 


( do 




-B-zw + D = 0 


(5.4) 


As mentioned, we will retain the DHD approximation (eq. (2.5)) after the stall, and we will continue to 
plot the ratio f? DHD approx °f neglected to retained terms by using the formula 


„ stalled PIO 
"DI ID approx = 


1 


d d V 

D 

B- 

— + v,. — In a 
\dt dr) 


(5.5) 


The same idealized inner velocity fields (eqs. (3.1)) are used in the analysis after the stall, and the 
momentum equation (3.3) can be written 


-2 Vp = pA 


(5.6) 
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where 


A] = Gi — G 3 V + Du 

(5.7a) 

A 2 = G'-> + Gt,u + Dv 

(5.7b) 

2 ■ 

A 3 = 2g + — D - 4Q y u 
e 

(5.7c) 

and G is given by equations (3.5). Substituting the density (eq. (5.3)) into equation (5.6), enforcing the 
cross-derivative equations (3.6), and defining 

L = In 0 - eZ 

(5.8) 

we obtain 


A x VL = curlA 

(5.9) 

By equation (5.6), we have 


£ 

X 

ClL3 

<N 

1 

II 

< 

Th 

3 

O 

(5.10) 

Although the flow at the interface r = a is baroclinic because of the jump in density, we take the inner 
flow to be barotropic so that 

curlA = 0 (r < a) 

(5.11) 

and 


A x VL = 0 (r < a) 

(5.12) 

If equations (5.7), (3.5), and (3.1) are substituted into equation (5.11), we obtain 


%= l -[^ + AQ z )-Dl]-n y l 

(5.13a) 

ri = -^[|(t + 4Qj + Lhi]-Q y L> 

(5.13b) 

t = -D(t + 2Q Z ) 

(5.13c) 


The expression (eq. (5.13c)) for t, is identical to equation (3.7c) for the PIO. Except for the Q v terms, 
equations (5.13a) and (5.13b) for f and f] are the same as equations (3.7a) and (3.7b) for the PIO. 
However, after a stall, there are no equations for D, x c , or y c that are comparable to equations (3.8), 
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(3.9a), and (3.9b) for the PIO. The absence of the D equation is fortunate because we now need D to be 
determined solely by the prescribed buoyancy B(t). The absence of equations for x c and y c means that 
the inner pressure can be determined for any given values of x c and y c ; that is, the track after a stall is 
not determined by our midtropospheric model as it was before the stall had occurred. 


When results (eqs. (5.13)) are substituted back into equation (5.7), we obtain 


2D 


A x =aX -2Q y D Z D (x c -U h ) + t(y c -V h ) + r\-4Q z V b 


(5.14a) 


where 


2D. 


A 2 =aY + 2Q y XZ-t,{x c -U b )-D(y c -V b ) 1 + 4Q z U b 


At, = 2 


[ D 

- 4Qy 

8 + — 

\ e ) 



U h+r[ Z + ^(DX-^Y) 


a (t) = D~ — 
2 


^ + 4Q Z )-D 2 


(5.14b) 

(5.14c) 


(5.15) 


Note that a(t ) is a measure of the departure from the PIO after the stall because if equation (3.8) for the 
PIO were enforced, a(t) would vanish. 

Equation (5.12) is satisfied by taking 


VL = C x {t) A 


(5.16) 


Integrating this equation for L gives 


L = Ci 


-( X 2 +Y 2 ) + 
2 ) 


-D{x c -U b ) + t,(y c -V b ) + ——r\-4 £i z V b + 4Q V — 




2D, 


-Z(x c -U b )-D(y c -V b ) — % + 4Q z U b 


Y + 2\g + --2Q U b \Z 


+ 2Q, y (-DX + tY -i\Z)Z 


+ 


(5.17) 


According to equations (5.3) and (5.8), the inner density can be written 


P = (!-£) e 


(5.18) 
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or 


p = C 6 (?)(l -B)exp 


Ci 


f(* 2 ^) + 


-D(*c ~ Ul , ) + 1 ( % - V b ) + — r, - 4Q .V b + 4Q 

e e 


2D, 


-Z(x c -U b )-D(y c -V b ) 1 + 4 Q z U b 


+2Q. y {-DX + XJ -r\Z)Z 


Y + 2 


D 

8 + 2Q U b 

e 


(5.19) 


The inner pressure p can now be obtained by substituting equations (5.19) and (5.14) into equation (5.6) 
and integrating to obtain the barotropic result 


1 

P 

2 Cj(0 


(5.20) 


where p is given by equation (5.19). This intermediate result confirms that equations (5.13) are 
acceptable ODEs for the inner vorticity of a stalled-PIO model that can translate with any given velocity 
[x c (r),y c (t)]- The next step is to obtain the outer solution and apply the jump conditions at the interface 

to obtain equations for a and D. As with the PIO model, we revert to a simplified axisymmetric, 
nontranslating solution that is sufficient for obtaining these two equations. 

6. Axisymmetric, Nontranslating Stalled-PIO Solution 

6.1. Inner Solution 


The inner density (eq. (5.19)) and pressure (eq. (5.20)) for the stalled PIO become axisymmetric and 
nontranslating by setting x c = y c =U b = V b =’E, = r\ = Q y =0. The inner velocity fields in cylindrical 

coordinates are again given by equations (4.1), and the inner density becomes 


p(r,z,t) = C 6 (t)(l-B) expjC](t) 


—ar 2 + 2 
2 


8 + ~ ■ 
V e 


Applying the boundary condition (eq. (5.2)) at r = a gives 


Ci(0 = - 


/ 

<? + — 

\ e / 


( 6 . 1 ) 


( 6 . 2 ) 


C 6 (t) = a b (t) exp 


eaa 


2 


8 + 


\ 


D N 

e / 


(6.3) 
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and the inner density (eq. (6.1)) and pressure (eq. (5.20)) become 


p(r,z,t) = a^(0(l-S) exp 


-eZ + - 


1 + * 

eg 


( r \ 

2' 

1- - 


\a) 



p(r,z,t) = -o®(t)(\-B) 


+ 

l £ SJ 


exp 


-eZ + - 


D 

eg 


( r \ 

2' 

1- - 


\a) 



where (5 is dimensionless and defined by 

2 


2 r 


n/ . £ a . . ea 

P (0 = — «(0 = — 

4 g 8g 


2Z)-^(C + 4Q-.) + D 2 


(6.4) 


(6.5) 


( 6 . 6 ) 


6.2. Outer Solution 

The outer velocity field for the stalled-PIO model is again given by equations (4.3) and the equation 
for a by equation (4.4). An equation for D can be determined from the jump condition that the pressure 
be continuous at r = a. To determine the outer pressure field, we solve the momentum equation (3.3) for 
the pressure gradient in cylindrical coordinates and then substitute the outer velocity field from equa- 
tions (4.3), t from equations (5.13c), and a from equation (4.4) to obtain 


. 0 

dp _ gp 

f( a ) 
y-F - 

2' 

dr er 

\rj 



(6.7a) 


dp 


0 


dp 

dz 


ae 

0 


= 0 


: -gp 


(6.7b) 

(6.7c) 


where F is given by (4.6) and 


y(t) = ^-(D-2Q£+D 2 
2 § y 


( 6 . 8 ) 


For the PIO, F was a constant of the motion; however, after the stall, F becomes time-dependent. 
Again we impose equation (4.7), so the outer density p° must satisfy 


dp u 

dr 


er 


y-F 


dp 

dz 


0 


(6.9) 
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The solution of this equation that satisfies the boundary condition (eq. (5.1)) at r = a is given by 


p °(r,z,t) = o b 


( °v 

„ F 

- exp 

-£Z + — 


2 



(6.10) 


where we must have y>0 for p°(r,z,t) to remain bounded as r — > oo. Substituting this result in 
equation (6.7) and integrating for the outer pressure, we obtain the hydrostatic result 


P°(r,z,t ) 


-p°(r,z,f) = o°-[- N J exp< 
e ^\r 


-eZ + - 


1- 



(6.11) 


For the PIO, p°(r,z,t) and p°(r,z,t) had maximum values at r = a\ however, for the stalled PIO, 
these maxima occur at radius r max , where 


f F \ 1/2 


\y ) 


(6.12) 


Again we identify these maximum values P° (r max ,b} and p 0 (q na x’^) at height b with the environ- 
mental pressure pf"'"' and environmental density p/) nvu at height b, which are constants. If we sub- 
stitute the radius given by equation (6.12) into the outer pressure (eq. (6. 1 1)) at height z = b, we can solve 
for o b (t) 


f F' 

\2 

exp 

i(y-F) 

II 

sJ-O) 

3 

< 

1 " 

'f' 

2 

exp 

I(y-f) 

\ Y 


L2 

8 

,y ) 


[2 


a^) = pr ir 

and the outer pressure (eq. (6.1 1)) and inner pressure (eq. (6.5)) become 


p"(r,z,t) = ±p' J (r,z,t) = p c b " 


'f' 

>~\ (N 

„ 1 

Ja \ 2 1 

^y) 

\7j exp 

-eZ + — 
2 

r W 1 


and 


p = pt wl \\-B) 


fl + ^l 


A, 


Y 


exp 


-eZ + — (y - F) + — ^-r- 
2 , D 

1 + — 

zg 


i- 

-fl 


\a) 


(6.13) 


(6.14) 


(6.15) 


6.3. Buoyancy and D 

By equating the outer pressure (eq. (6.14)) and the inner pressure (eq. (6.15)) at r = a, we again obtain 
the result (eq. (4.15)); but now, after the stall, the buoyancy B(t) is given, and w (or D) is determined 
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by integrating this equation. After the stall, equations (4.22) and (4.23) for the outer and inner tempera- 
tures also remain valid. With equation (4.15) substituted, the inner pressure (eq. (6.15)) becomes 


p(r,z,t) = p| nvir — I exp J -eZ + — (y - F) + A. l-(- 


(6.16) 


2 r 

Ht) = (3(1 - B) = i-\2egB - (1 - B) £(£ + 4Q Z ) - D 2 
8s l 4 


(6.17) 


and now 


2 

y(t) = ^-UgB-2QX + D : 
2p \ 


(6.18) 


We substitute equations (6.16) and (4.23) into equation (4.25) and integrate over the volume of the core to 
obtain the thermal input power 


-v stalled PIO envir ^ (X|2 | [l 


exp — (y -F) Nl-exp[-e(/i-&)]} 




? d: I 2\F ) 2 \y ) X X : 


e x (l-X)-l l (6.19) 


By equations (6.4), (4.15), and (6.17), the inner density p is given by 


Y r 

p(r,z,t) = -Q.-B)p(r,z,t) = p/® nvir (l-5)-[ — 2 exp- -eZ + -(y-F) + X l-(- 

8 8\ y 2 \a 


( 6 . 20 ) 


and by comparison with equation (5.3) 


o(r,t) = P/f vir - f— f exp U 




( 6 . 21 ) 


Substituting equation (6.21) into equation (5.5) gives for the ratio 


D stalled PIO , , 1 \ u F Y , 1 Y , ( F < , [ r 

« D HD,pp ra (>-.0-- - - 


( 6 . 22 ) 
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Because this ratio depends upon r, we will evaluate it during model runs at r = 0, r = 0.707a, and r = a 
and choose the maximum (i.e., conservative) value. The cloud-base mass influx after the stall is obtained 
by substituting equation (6.20) into equation (4.19). 


^stalled PIO 


-a 2 D(\-B)pl WV[ \ 

'f s 

|2 

exp 

» 1 ^ 
i 

i 

8 

\ Y ) 


L2 J 


(6.23) 


6,4. Central Pressure Deficit 

The environmental pressure p envir (z) is obtained as a function of height by substituting equa- 
tion (6.12) into equation (6.14) 


envir/ \ envir -eZ 

P ( z) = Pb e (6.24) 

We find the minimum internal pressure p m [ n (z,t) at the center of the mesocyclone, which is now 
becoming a tornadic cyclone, by setting r = 0 in equation (6. 16) 


p mm (z,t)= P e r r 


Y 


(f\ 2 

, 1 

— exp 

-eZ + A + -(y-F) 

\Y / 

2 ' 


The central pressure deficit A p(z,t) is defined as 


or 


A p(z, t ) = p envir (z) - p min (z, t) 


Ap(z,t) = pf nvir e~ eZ I 1 


( F ) 

2 

ft 

- 

exp 

A + -(y-F) 

UJ 


2 


(6.25) 


(6.26) 


(6.27) 


where A, y, and F are given by equations (6.17), (6.18), and (4.6), respectively. In our calculations, we 
will always evaluate A p at height b where Z = 0. 

7. Tornadic Cyclone Solution 

7.1. Constant B After Stall 

Figure 4 shows the stalled-PIO solution for Tornadic Cyclone A9, where the truncated value of B is 
held constant after the stall that starts at t =* 21 min. The environmental parameters are given in table 1 
and the initial values in table 3. (Purely for pedagogy, this figure also shows the plots for a fictitious 
unstalled PIO where A T unrealistically reaches +96 K.) The truncated buoyancy shown in figure 4(a) is 
insufficient to cause the rapid draft reversal required by the PIO (fig. 4(b)), so the contracting downdraft 
is prolonged, as shown in figures 4(b) and (c). Prolongation of the contracting downdraft spins up the 
Rossby number well above the tornadic cyclone value of Ro « 1000, as shown in figure 4(d). The ratio 
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^DHD approx from equations (2.7) and (6.22) is plotted in figure 4(e). This ratio is small, except for a 

brief interval when it reaches 0.28 just before the stall begins. Therefore, the DHD approximation was 
reasonably accurate during this segment of the calculation. 

The prolonged downdraft of the stalled-PIO model is in qualitative agreement with the collapsing 
phase of a supercell storm, as described by Rotunno (1986): “The BWER [bounded weak echo region] 
begins to fill, and downdrafts intensify ... . In association with this development, a tornado forms at full 
strength and may last from a few to several tens of minutes.” The prolonged downdraft is also in agree- 
ment with more recent observations by Trapp (2000) of weak-to-moderate downflow throughout the 
entire depth of a mesocyclone just before tornadogenesis. It is also in agreement with early observations 
by Fujita (1972) that tornadoes are most likely to occur during a decrease in the height of clouds that 
overshoot the top of a supercell. 

The mesocyclonic circulation T is plotted in figure 4(f), where it again achieves the value 
5 2—1 

r ma x =5x10' m“ s . The initial values for C, and D in table 3 were chosen so that T = 0 at t = 0. 
The cloud base mass influx M is plotted in figure 4(g) and the thermal input power Q in figure 4(h). 
Both of these quantities go through zero at maximum spin-up when w goes through zero. 

The central pressure deficit is computed from equation (6.27) at height z = b and is plotted in fig- 
ure 4(i). Before the stall, we see that A p = 0, which is consistent with inertial flow. After the stall, A p 
builds to a value that substantially exceeds the maximum measured value of 34 mb. The radial depend- 
ence of the pressure at height z = b is plotted in figure 4(j) before the stall and at various times after the 
stall. These plots were computed from equation (6. 16) for r < a and from equation (6.14) for r > a. 

1 . 2 . Extending Lifetime of Simulated Tornadic Cyclone 

An obvious shortcoming of our simulated tornadic cyclone is that it lasts for only about 4 min, while 
Tornadic Cyclone A9 lasted for 80 min. According to figures 4(b)-(d), the simulated tornadic cyclone 
began to expand and spin down when w became positive at t ~ 25.7 min. However, should we prescribe 
that the buoyancy B -* 0 as w —> 0, then by equations (2.5), (4.15), and (5.13c), it follows that 
D, w, D, and t, —» 0 also, and the simulation would stay in a state of maximum spin-up indefinitely. 
(Recall that for an actual storm, B = w = 0 means that the average buoyancy and average vertical flow 
would remain zero; that is, the updrafts and downdrafts would remain in balance.) 

We can try to apply this intuitive argument to our mathematical solution. However, should we set 
B = w = D = D = t, = 0, we immediately encounter two problems: (a) as given by equation (6.18), y 

would then be negative and our solution would lose its validity because the outer pressure (eq. (6.14)) 
would become unbounded as r -*■ oo; and (b) the DHD approximation would become inaccurate because 

the ratio 7 ?dhd approx’ as given by equation (6.22), could become very large. To solve problem (a), we 
will assume that when Ro > 1000, we may neglect the Coriolis force by letting Q,. — >0 so that the trou- 
blesome term -2Q X, disappears from equation (6.18). The transition from inertial to cyclostrophic flow 
that occurs after the stall supports this assumption. We can ameliorate problem (b) by specifying that as 
B -» 0, D -* -1.875 x 10 -5 s _1 instead of going to zero, which corresponds to w —* - 0.158 m s -1 . Al- 
though the core continues to contract, Ro increases by less than 10 percent during the 80-min simulated 
lifetime, and the DHD approximation remains reasonably accurate. 
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Implementing this approach after the stall, we replace the constant B with 


and the constant Q z with 

Q' z (t) = Q z [l-H n {t-t 0 )\ 

where H n (t) is a smoothed Heaviside unit step function defined by 

H n {t) = ^[l + tanh(nt)] 


(7.1) 


(7.2) 


(7.3) 


The values n = 0.011111 and f 0 =1543.7 s provide a smooth transition and achieve the target value 
for D (or w) given previously. The resultant plots for a sustained tornadic cyclone are shown in 
figure 5. During the sustained period (t>30min), the values for M and Q are relatively small 

( M = -6x10 kg s and Q = -3 x 10 W) because the updrafts and downdrafts are nearly in balance. 


7.3. Terminating the Tornadic Cyclone 

We shall now attempt to simulate the observed termination of Tornadic Cyclone A9 after 80 min. Our 
strategy is to use B(t ) to induce a positive pulse in D(t) or wit) that will cause the core to expand and 
spin down, after which w also returns to zero. This approach requires Bit) to make a positive excursion 
followed by an equal negative excursion. But again we encounter problems: (a) because y goes negative 
during the negative excursion of B(t) and (b) because approx has an infinite pole where D(t) goes 

through zero at the onset of its positive pulse. 

We can prevent the problem with y by decreasing the negative excursion of Bit), which has the con- 
sequence of leaving D (or w) in a somewhat elevated or nonzero final state — which turns out to be for- 
tuitous. Because Rdhd approx* as given by equation (6.22), is a function of both r and t, we cannot 
circumvent its pole at t = 103 min by further tailoring B(t). However, we can argue that the ratio 
T^DHD approx as defined by equation (6.22) could be too severe a test of the DHD approximation, because 

anytime that D passes through zero, we risk such a pole. Regardless, we will proceed so that the model 
will conform to the data. 


After the stall, the complete buoyancy is given by 


B"(t) = 


12 


288.6 


1 8 

( — (t — h) ) + 


H m (t- h )-3Hi(t-t l ) + 2H 3 m (t-t l )] [l-ktf,(t-T)] 


(7.4) 


where n = 0.011111, m = 0.005, s = 1, /p = 1543.7 s, /y = 6761 s, and k = 0.7432. On the right-hand 
side, the first two terms are the same as in equation (7.1), the third degree polynomial in H m (/ - ?| ) gives 

the terminating bipolar excursion, and the right-most factor decreases the negative part of this excursion. 
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For further discussion of polynomials and generalized functions based on the hyperbolic tangent form of 
H n (t), see Costen (1967). 

The complete simulation of Tomadic Cyclone A9 is plotted in figure 6. It is apparent from figure 6(a) 
that we have sought a small perturbation in buoyancy that would be effective in terminating the tornadic 
cyclone. An enlargement of this perturbation is shown in figure 7(a). The effects on w and a are shown 
in figures 6(b) and 6(c). The final elevated value of w causes an increase in a and dramatic declines in 
both Ro and A p, as shown in figures 6(d) and 6(e). The ratio R^HD approx ' s plotted in figure 6(f). The 

pole at t » 103 min coincides with the final ascent of w through zero. The discontinuity in approx 

at t « 1 1 3 min results from the discontinuity in B that occurs between the positive and negative lobes of 
A T, as shown in figure 7(a). 

The mesocyclonic circulation T, as shown in figure 6(g), is unaffected during termination, which con- 
firms that the declines in Ro and A p are the result of core expansion from the induced updraft. The 
mass influx at cloud base M is plotted in figure 6(h). The increase in M that results from the expanding 
updraft is seen by enlarging the last 40 min, as shown in figure 7(b). Near the end, the computed value 
for M is rapidly increasing and, if the computation continued further, M would soon reach the mea- 
sured reference value M re f = 1.2 x 10 9 kg s _1 reported by Foote and Fankhauser (1973) for a different 
supercell. The thermal input power Q is shown in figure 6(i), and a similar enlargement is shown in 
figure 7(c), where the discontinuity in Q again marks the end of the positive lobe of AT. Like M, Q is 
rapidly increasing, and if the computation continued further, Q would soon reach the measured reference 

value Q ie f = 1.9xl0 13 W. 


During the prolonged phase of the model tornadic cyclone, AT = 0 and w = -0.158 ms 1 . Since 
these quantities represent average values for an actual storm, we infer that storm longevity requires the 
updrafts and downdrafts to be essentially equal in magnitude and the positive buoyancy that drives the 
updrafts to be in balance with the negative buoyancy that drives the downdrafts. According to the model, 
termination of the tornadic cyclone requires a thermal input that warms and intensifies the updrafts, 
warms and diminishes the downdrafts, or both. In nature, this thermal input would result from increased 
condensational heating, decreased evaporative cooling, or both. This model-based concept of weakening 
a tornadic cyclone by the application of heat is consistent with the observed weakening of tornadic storms 
that move from land to water, which is a source of additional heat. 

Following Kessler (1972) and others, we speculate how human beings might intervene to trigger the 
termination process earlier. According to figures 7(a) and 7(c), the positive lobe of AT results from a 
thermal power input of 11 -min duration that peaks at 158 GW. The average power input is 61 GW, 
and the total thermal energy deposited into the tornadic cyclone is 40 TJ. As a result of this thermal input, 
w increases from -0.158 m s-1 to +10 m s -1 . The negative lobe of AT subsequently reduces w to 
+8 m s— !, and this sustained average updraft causes the core to expand and spin down quickly. 

Suppose that we were able to manually input 40 TJ of thermal energy into the tornadic cyclone at 
some earlier time with the result that w is increased to 10 m s-f In a convectively unstable troposphere, 
this new value for w could be sustained naturally by latent energy release without further manual input. 
If the new value for w were sustained, we would have successfully triggered the termination process. 
Selectively injecting this thermal energy into the downdrafts would avoid augmenting the production of 
hail. 


24 



Our value 40 TJ (or 11 million kWh) is an order of magnitude lower than Kessler’s (1972) estimated 
input energy of 80 million kWh to prevent the initial development of a supercell storm. At a rate of 
0.1 USD per kWh, 40 TJ would cost 1 million USD. This expenditure compares favorably with the 
1 billion USD property damage plus injuries and loss of life incurred by Tornado A9. 

The model also provides insight into forecasting which supercell storms will develop a tornadic 
cyclone and produce strong tornadoes. If the supercell conforms to the PIO model, it can produce hail but 
not strong tornadoes. Therefore, if Doppler radar images show that parcels in the core are following iner- 
tial trajectories (as projected on a horizontal plane), the supercell will not produce a strong tornado. 
However, should the core trajectories transition from inertial to cyclostrophic, the forecaster can infer that 
the PIO has stalled and that the formation of a tornadic cyclone and a strong tornado is likely. 

8. Concluding Remarks 

Unlike the PIO model or the conventional splitting-storm model, the stalled-PIO model is capable of 
simulating a tornadic cyclone that originates in the midtropospheric layer — in agreement with radar ob- 
servations of the strongest tornadic storms. After the PIO model became stalled in a state of contraction 
and spin-up, we were able to tailor the buoyancy in the core so that the stalled-PIO model simulated the 
80-min lifetime of intense Tornadic Cyclone A9 that occurred over Oklahoma City, Oklahoma, on May 3, 
1999. 

The dilatation-horizontal divergence (DHD) approximation used in the PIO and stalled-PIO models 
states that the horizontal divergence in the cylindrical core is predominantly due to the three-dimensional 
dilatation of rising/falling air parcels. Thus, the horizontal divergence is approximately equal to the up- 
draft speed divided by the scale height. To check the accuracy of this approximation during our model 
runs, we plotted the ratio of neglected terms to retained terms in the continuity equation. This check was 
quite severe because whenever the horizontal divergence passed through zero, the plotted ratio could 
blow up. Such a blowup occurred once, but it was near the end of the run when the tornadic cyclone was 
decaying. Twice the ratio approached the value 0.3, although most of the time the ratio was much less 
than 0.1. We conclude that the DHD approximation was reasonably accurate during our PIO and stalled- 
PIO model runs. 

The model tornadic cyclone remained in a nearly steady state when the spatial average buoyancy, 
average vertical flow, and average thermal input from latent energy release were close to zero, that is, 
when the updrafts and downdrafts were in balance. The tornadic cyclone was terminated by a predomi- 
nantly positive pulse of buoyancy that increased the average vertical speed to 8 m s-f As a consequence 
of the DHD approximation, this sustained updraft caused the tornadic cyclone to expand and spin down 
rapidly. 

The terminating buoyant pulse required a thermal input of 40 TJ or 1 1 million kWh. For Tornadic 
Cyclone A9, this thermal input was supplied by latent energy release. However, we could, in principle, 
manually inject this thermal energy to trigger an early termination of the tornadic cyclone and its pendant 
tornado. Our calculated energy value of 40 TJ is an order of magnitude less than previous model-based 
estimates for effective human intervention. 
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Table 1. Parameters Used in PIO and Stalled-PIO Models for Tornadic Cyclone A9, Where Meteoro- 
logical Values Are From Edwards and Thompson (2000) and Burgess and Magsig (2000). 


Qj., rad s _1 
b, km MSL 
h, km MSL 

Pb > kra 


T°(b), K 

dT°/dz, K m _1 
g, m s -2 
R, Jkg -1 K _1 
c p , Jkg -1 K- 1 
L c , J kg -1 . . . 


4.114 x 10 -5 

1.727 

7.810 

82 

1.1842xl0 -4 

288.6 

.. -8 x 10 -3 

9.81 

287 

1004 

. . . 2.5 xlO 6 


Table 2. Initial Values Used in Maximal PIO Model Run (Figs. 1 and 3), Where (-12 K < AT < 12 K) 

and r max = 5 x 10 5 m 2 s _1 


£(0), s" 1 -5.0748 x 10 -5 

D(0), s" 1 -5.4972 xl0~ 4 

a( 0), km 71.344656 


Table 3. Initial Values Used in Fictitious PIO Model (Fig. 4), Where (-12 K < AT <96 K), and Stalled- 
PIO Model (Figs. 4 to 7), Where (-12 K < AT < 12 K) and r max = 5 x 10 5 m 2 s -1 


^(0), s" 1 -1.03956 xl0~ 8 

0(0), s _1 -1.50804 x 10 -3 

a( 0), km 43.9939 
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(a) Contraction and spin-up during first hour. 



(b) Two-hr trajectories. 

Figure 1. Aerial view showing inertial trajectories of parcels on core periphery during 2-hr period centered on 
maximal PIO pulse. (See also fig. 3.) 
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Figure 2. Pulsing inertial oscillation (PIO) model or stalled-PIO model for convective core of supercell storm or 
tomadic cyclone. The core radius is a(t) and centerline is located at [x c The horizontal divergence 
D(t), buoyancy B(t), and vertical components of vorticity E,(f) and velocity w(t) are uniform inside core and 
zero in outer region. The midtropospheric layer extends from heights b to h. 
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(a) Core temperature excess AT versus t. 



(b) Updraft speed w versus t. 

Figure 3. Maximal PIO plots during 2-hr period centered on pulse. (See also fig. 1.) As in all runs, mesocyclonic 

5 2—1 

circulation at maximum contraction and spin-up has value r max = 5 x 10 m ~ s . 
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(c) Core radius a versus t. 



(d) Rossby number Ro versus t. (See eq. (3.2).) 
Figure 3. Continued. 
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(e) Ratio ^dhd approx °f neglected terms to retained terms in dilatation-horizontal divergence (DHD) approximation 
(eq. (2.5)) versus t. (See eq. (2.7).) 


x 10 5 



(f) Mesocyclonic circulation T versus?. (See eq. (4.18).) 
Figure 3. Continued. 
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(g) Cloud base mass influx M versus t. (See eq. (4.20).) 

x io 14 



(h) Thermal input power Q versus t. (See eq. (4.26).) 
Figure 3. Concluded. 
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(a) Core temperature excess AT versus t. 



(b) Updraft speed w versus t. 

Figure 4. Stalled-PIO plots for Tornadic Cyclone A9, with buoyancy B held constant after stall. Also shown for 
pedagogy are plots for fictitious unstalled PIO. 
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(c) Core radius a versus t. 



t, min 

(d) Rossby number Ro versus t. (See eq. (3.2).) 
Figure 4. Continued. 
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(e) Ratio R DHD approx °f neglected terms to retained terms in DHD approximation (eq. (2.5)) versus t. (See 
eqs. (2.7) and (6.22).) 


r, m 2 s~ 


X 10 5 



(f) Mesocyclonic circulation r versus?. (See eq. (4.18).) 
Figure 4. Continued. 
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X 10 10 



(g) Cloud base mass influx M versus t. (See eqs. (4.20) and (6.23).) 
x io 14 



(h) Thermal input power Q versus t. (See eqs. (4.26) and (6.19).) 
Figure 4. Continued. 
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(i) Central pressure deficit A p versus t. (See eq. (6.27).) 



(j) Pressure p versus radius r before stall and at various times after stall. (See eqs. (6.14) and (6.16).) 

Figure 4. Concluded. 
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(a) Core temperature excess A T versus t. 



(b) Updraft speed w versus t. 

Figure 5. Stalled-PIO plots for sustaining Tornadic Cyclone A9 by decreasing the buoyancy B and gradually 
neglecting the Coriolis force after stall. 
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(c) Core radius a versus t. 



(d) Rossby number Ro versus t. 
Figure 5. Continued. 
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(e) Ratio R DHD approx °f neglected terms to retained terms in DHD approximation (eq. (2.5)) versus t. 
x 10 5 



(f) Mesocyclonic circulation T versus t. 
Figure 5. Continued. 
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(g) Cloud base mass influx M versus t. 



(h) Thermal input power Q versus t. 
Figure 5. Continued. 
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(i) Central pressure deficit A p versus t. 
Figure 5. Concluded. 
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(a) Core temperature excess A T versus t. 



(b) Updraft speed w versus t. 

Figure 6. Final stalled-PIO plots for Tornadic Cyclone A9, including its termination by further tailoring of 
buoyancy B(t). 
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(c) Core radius a versus t. 



(d) Rossby number Ro versus t. 
Figure 6. Continued. 
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(e) Central pressure deficit A p versus t. 



(f) Ratio Rqhd approx °f neglected terms to retained terms in DHD approximation (eq. (2.5)) versus t. 

Figure 6. Continued. 
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(g) Mesocyclonic circulation T versus t. 



(h) Cloud base mass influx M versus t. 
Figure 6. Continued. 
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(i) Thermal input power Q versus t. 

Figure 6. Concluded. 
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(a) Core temperature excess A T versus t. 
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(b) Cloud base mass influx M versus t. 

Figure 7. Enlarged plots of final stalled-PIO during termination phase of Tomadic Cyclone A9. 
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(c) Thermal input power Q versus t. 


Figure 7. Concluded. 
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